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Abstract. The column density distribution of trapped OH~ ions in a 22-pole 
ion trap is measured for different trap parameters. The density is obtained from 
position-dependent photodetachment rate measurements. Overall, agreement is 
found with the effective potential of an ideal 22-pole. However, in addition we 
observe 10 distinct minima in the trapping potential, which indicate a breaking of 
the 22-fold symmetry. Numerical simulations show that a displacement of a subset 
of the radiofrequency electrodes can serve as an explanation for this symmetry 
breaking. 
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1. Introduction 

Multipole radiofrequency ion traps [1], in particular the 22-pole ion trap [2], are 
versatile devices used in laser spectroscopy [3j HJ [5l El El O [9l [10] and investigations of 
chemical reaction processes [H [11] [T^ [T3] of atomic and molecular ions. High order 
multipole traps offer a large field free region in the trap center, and therefore provide 
a reduced interaction time of the ions with the oscillating electric field compared 
to a quadrupole trap [14] . Buffer gas cooling down to cryogenic temperatures of 
the translational [31 [7] , rotational [21 [S] , and vibrational degrees of freedom of 
trapped molecular ions has been demonstrated. This enables applications with all 
the advantages of low temperature experiments, such as a reduced Doppler width in 
spectroscopy studies and a well defined population of internal states. 

Stable confinement of a single ion in the oscillating quadrupole field of a Paul 
trap is precisely predicted, because the Mathieu equations of motion can be solved 
analytically. This is different for oscillating high order multipole fields, where the 
equations of motion have no analytical solution. However, the movement of the ions 
in a fast oscillating rf-field justifies the assumption of an effective trapping potential, 
based on the separation of the ion motion into a smooth drift and a rapid oscillation, 
called micromotion [15]. This effective potential can be expressed as 



with charge q in an electric field oscillating on frequency f2 and $o is a non-oscillating 
dc potential. For an ideal cylindrical multipole of the order n this effective potential 
can be expressed as 



where Vq denotes the amplitude of the oscillating rf-field. For a high order multipole 
field (n = 11 for the 22-pole trap) this creates an almost box- like trapping volume with 
steep walls and a large field free region in the center. This is important to suppress rf 
heating when cooling ions in neutral buffer gas [T^ . 

The number of trapped ions that can be excited with radiation depends on the 
local density of ions in the interaction volume with the light field. In experiments 
with trapped Ba"*" ions in an octupole trap the density has been imaged by spatially 
resolving the fluorescence signal [Ml [ITj. In this article we report on a method 
based on the photodetachment of trapped anions that allows us to image the density 
distribution of the trapped ions. In this way we measure the experimental trapping 
potential of a 22-pole trap with unprecedented accuracy and study its dependence 
on the radiofrequency field. In the next section the photodetachment tomography is 
explained, followed by a presentation of measured density distributions at different 
temperatures and rf amplitudes. We find ten unexpected maxima in the density 
distribution and analyze their origin with the help of numerical simulations in section 



where E{r) denotes the electric field at the point r, m is the mass of a test particle 





(2) 



2. Experimental procedure 



The experimental radiofrequency ion trap setup has been described in detail in Ref. 
[18] . Its central element is a 22-pole trap [2] with good optical access along the trap 
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axis. In the present tomography experiments, we have used an ensemble of stored 
0H~ anions cooled to a translational temperature of either 170 or 300 K. The use of 
negative ions allows us to perform photodetachment measurements, were the anions 
are depleted by a laser beam propagating parallel to the symmetry axis of the trap 
[191 [20] . The photodetachment process 

OH" + ht^ ^ OH + e" (3) 

with hv > 1.83 eV yields a depletion rate proportional to the photodetachment cross 
section and to the overlap of the ion column density with the photon flux. By scanning 
the position of the laser beam and measuring the depletion rate at each position the 
relative ion column density distribution is obtained. For small ion densities, where 
the Coulomb interaction between the ions can be neglected, the depletion rate at the 
transverse position {x, y) can be expressed as 

k{x, y) = (TpdFLp{x, y). (4) 

It depends only on the total photon flux Fl, the photodetachment cross section apd 
and the single particle column density p{x, y), reflecting the spatial overlap of the laser 
beam with the ion distribution I19j . As the latter is normalized to unity, we can write 
it as 

jgk(x,y)dS 

Thus, a full two-dimensional tomography scan of the trapped ions can be used to map 
the entire ion column density in the trap. If k{x,y) is measured for the whole ion 
distribution, also the absolute photodetachment cross section can be obtained 20j . 

Ions are produced in a pulsed supersonic expansion of a suitable precursor gas, 
crossed by a 1 keV electron beam. For the OH" production we use a mixture of 
Ar/NH3/H20 (88%, 10%, 2%) to create NH^. A rapid chemical conversion by water 
forms the OH" anions in the source. A bunch of ^--^500 mass-selected ions is loaded 
into the trap, which is enclosed by a copper housing that is temperature-variable 
between 8 and 300 K. A typical He buffer gas density of 1 x 10^'*cm~'^, employed 
for all temperatures, is enhanced by a buffer gas pulse during the ion injection. The 
trap is operated with fl = 2tt x 5 MHz radiofrequency and different rf amplitudes. 
Along the axial direction ions are confined by end electrodes biased to -2 V. To allow 
good thermalization of the ions, a storage period of 200 ms is inserted before the laser 
beam is switched on. After a given storage and laser interaction time the current 
signal, proportional to the number of ions that survived the interaction with the 
photodetachment laser, is detected with a microchannel plate. 

The two-dimensional tomography scans are performed under the same 
experimental conditions as in Ref. |20j . A free-running continuous wave diode laser 
at 661.9 nm (Mitsubishi ML1J27, lOOmW, spectral width 0.7 nm FWHM) has been 
employed. It is imaged into the trap and scanned in the (x, y) plane, perpendicular to 
the symmetry axis of the trap, by moving the imaging lense on a mesh with 0.25 mm 
point spacing with a computer-controlled two-dimensional translation stage. The mesh 
spacing is comparable to the 1/e^ radius of the laser in the trap of 350 ^m. For each 
laser position the photodetachment depletion rate is obtained in a storage time interval 
of 0.2-2 s by fitting an exponential decay to the ion current signal, reduced by typically 
one percent due to the background loss rate. The data are averaged over typically 
4-8 scans. Within each scan the mesh points are accessed in random order to avoid 
systematic drifts. 
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3. Tomography of the trapping potential 

Fig. [T^ shows a tomography scan of 0H~ anions in the 22-pole trap at 300 K with the 
rf amplitude set to 160 V. The figure also contains the sketched arrangement of the 
trap's copper housing mounted on the coldhead of the cryostat, the position of the 22 
rf electrodes and the axial end electrodes. Fig. [T|3 is a zoom of the scan which more 
clearly shows the measured ion density distribution. Every pixel of the histogram 
here represents a fitted photodetachment depletion rate k{x,y) (see previous section) 
and is proportional to the single-particle probability density p{x, y) along a column 
parallel to the z-axis. As can be seen the ion distribution as a first approximation 
can be considered circularly symmetric and constant in the center region of the trap, 
whereas it drops to zero when the ions reach the outer regions of the trapping volume. 
This distribution directly visualizes the overall storage properties of a 22-pole ion trap 
with a flat potential in the center and steep walls. Note that the ion density drops to 
zero already for smaller radii than the end electrode (solid line in Fig. [Hd), indicating 
that clipping of the laser at the end electrode is not affecting the measured density 
distribution. For smaller rf amplitudes the ion density distribution would extend to 
larger radii and could not be fully probed by the photodetachment tomography. For 
this reason we have restricted ourselves to large enough rf amplitudes in this study. 

In Fig. [TJ; a horizontal cut through the ion distribution is shown. While in the 
center the distribution is relatively uniform, the population is locally enhanced by 
up to 40% near the edge of the ion distribution. Such a behavior has already been 
observed in previous measurements [19j . To study this here in more detail, the effective 
potential V{x, y) is extracted from the local ion density p{x^ y) assuming a Boltzmann 
distribution for the ions in the trap 

P{xi y) ^ ^ exp{-V{x, y)/kBT), (6) 

where T is the absolute temperature, is Boltzmann's constant, and Z is the 
partition function. Since only the column density is measured, the resulting potential 
V{x, y) is an average over the z-dircction. Fig. [TJi shows a cut through the obtained 
effective potential for the distribution of Fig. [TJ:. Overall this potential compares well 
with the calculated potential of an ideal 22-pole potential (solid line) , obtained without 
any free parameters from Eq. Closer inspection reveals interesting features in the 
potential that deviate from the ideal multipole. While the potential is relatively flat 
in the center, it shows a distinct minimum of about 12meV near the left edge of the 
ion distribution and a weaker minimum of about 5 meV near its right edge. It can be 
excluded that this change of the distribution is caused by space charge effects, because 
the experiments are performed with only a few hundred ions in a trap volume of about 
1 cm^. 

The same features of the effective potential are observed in measurements at a 
lower trap temperature. Fig. [5^ shows a tomography scan at 170 K and the same rf 
amplitude as above. The ion distribution again looks circular symmetric with a distinct 
cutoff when the ions reach the steep walls of the trapping potential. A horizontal cut 
through the effective potential, obtained in the same fashion as Fig. [TJi, is shown in 
Fig. The same minima as for 300 K are observed in the effective potential. At this 
lower temperature the two minima are better resolved and appear similar in depth on 
the left side and slightly deeper on the right side of the potential as compared to the 
300 K tomography. 
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Further substructure becomes visible in the 170 K density distribution. Ten 
clearly separated maxima in the density distribution appear almost equally spaced 
in angle at a radial position of about 3 mm. According to Eq. ^ they correspond 
to ten localized minima in the trapping potential at this radius with a typical depth 
of lOmeV. These minima have not been significant in the 300 K ion distribution at 
160 V rf amplitude. They become visible, however, for larger amplitudes. Fig. [3^ 
shows a 300 K ion distribution for an rf amplitude of 270 V. It reveals the same ten 
density maxima and respective potential minima that could only be resolved at lower 
temperature at 160 V. 

We have studied the dependence of the depth of the ten potential minima on the 
rf amplitude at 300 K. The depth of the deepest minimum is plotted in Fig. [3}3. Since 
the effective potential is expected to depend quadratically on the rf amplitude, a fit 
with only a constant and a quadratic term is applied to the data (solid line in Fig. [Sh) . 
It yields an rf- independent offset of about llmeV, which is attributed to the static 
potential of the end electrodes of the ion trap. These end electrodes produce a radially 
repulsive potential inside the trap, as discussed in Ref. [191 HI] ■ It compares well with 
simulations, as shown in the next section. The ten "pockets" in the potential, however, 
reveal a more complex deviation from the ideal multipole description of Eq. ([2]). An 
explanation for them will also be discussed in the next section. 

4. Modeling trapping potentials of realistic multipoles 

The effective potential of a 2n-pole has a 2n-fold rotational symmetry, when averaging 
over one radiofrequency period. The appearance of the ten observed potential minima 
is therefore a clear indication for a breaking of the ideal symmetry. To investigate this 
effect further, the effective potential of the employed 22-pole trap has been modeled 
using a numerical simulation package based on a fast multipole solver [22l [23] for the 
boundary element problem in combination with accurate field evaluation in free space. 
With this method the electric field E{r) can be calculated at any location inside the 
trap. It is converted into the effective trapping potential using Eq. ([1]). We have 
verified that the simulation of the trapping potential of an ideal 22-pole structure 
reproduces the effective potential of Eq. ([2]) on the numerical level of accuracy. 

Different assumptions have been tested as origin of the observed ten potential 
minima, such as the influence of the shape and position of the end electrodes and of 
the copper housing around the trap electrodes, without showing a measurable effect 
on the potential. This suggests that imperfections of the trap geometry itself may 
be responsible. To simulate these imperfections a breaking of the ideal symmetry is 
introduced by displacing one half of the 22 radiofrequency electrodes by a small angle 
(see inset in Fig. |4^). Such a small tilt of one set of electrodes against the others 
occurs to be the most likely displacement during the trap assembly. Upon tilting one 
set of electrodes by only a few tenths of a degree the calculated effective potential of a 
22-pole trap at 160 V rf amplitude and -2 V on the end electrodes immediately shows 
ten potential minima. 

In Fig. [4^ the dependence of the maximum pocket depth on the tilt angle, as 
obtained from a series of simulations, is plotted. These simulations have been carried 
out for 160 V rf amplitude. Here, an imperfection in the parallelity of only 0.2° causes a 
pocket depth of 5 meV. The pocket depth is calculated at each angle with and without 
a potential of -2 V applied to the static end electrodes. The end electrode voltage 
produces an overall quadrupole potential that pushes the ion ensemble towards larger 
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radii in the trap in addition to the tih-induced pockets. Both data sets with and 
without end electrode potential are described by the same quadratic increase. For the 
simulations with end electrode potential a constant offset of about 9 meV is obtained, 
in fair agreement with the experiment value of about llmeV. 

From the measured depths of the potential minima (Fig. [8)3) a value of between 
3 and 5meV is extracted for an rf amplitude of about 160 V, after subtracting the 
influence of the static end electrode (see Fig. [Sb). Such an rf-induced pocket depth 
is obtained in the simulation for a tilt angle of between 0.15° and 0.2° (see Fig. H^). 
Fig. shows a simulated density distribution for a tilt angle of 0.15°. Following 
Eq. ([6]), the simulation has been performed for 0H~ ions that are stored at 300 K in 
the 22-pole trap with 270 V rf amplitude and -2V potential on the end electrodes. 
This simulated density distribution agrees well with the measured distribution of Fig. 
[3^, which has been obtained with the same trap parameters. A larger tilt angle was 
found to already overestimate the ten potential minima. Note that it is preferable to 
compare graphs of the density distributions of simulation and experiment instead of 
effective potentials, because the experimental potential is obtained by a logarithm of 
the density distribution which suppresses the fine details in the images. 

When the 22-pole trap was assembled the strong influence of small displacements 
of the rf electrodes on the effective potential was not known. A tilt of one set of rf 
electrodes by a tenth of a degree can therefore not be excluded for our trap. Such 
small tilt angles already come close to the mechanical tolerances for the assembly of a 
22-pole trap in the presently used design. This shows the need to significantly improve 
the precision in the rf electrode geometry when a potential with pocket depths in the 
/L^eV range is desired. 

In searching for an explanation for the observed ten minima, we have extended 
the electric field simulations to multipoles of different order n. These calculations have 
shown that the number of minima observed in the effective potential of a distorted 
multipole ion trap is directly connected to the multipole order as iVminima — n — 1. 
Besides the above discussed tilt of one set of electrodes, other distortions of the ion 
trap, such as a parallel displacement of one set of rf electrodes, also introduces n — 1 
minima. The number of minima therefore seems to be a general consequence of 
breaking the symmetry of a multipole ion trap. We expect the minima to be related 
to the points in space where the superimposed time-dependent electric field of the 
individual multipole electrodes cancels, because according to Eq. ([T]) these are the 
global minima of the effective potential. For a perfectly symmetric multipole trap 
cancellation is expected only in the trap center, but for a distorted symmetry several 
such points are found at larger radii. 

5. Conclusion 

In this article we report on a method to directly measure the column density 
distribution of ions in a 22-pole ion trap using photodetachment of stored OH". 
The two-dimensional tomography scans yield the effective potential averaged over the 
length of the trap. The measurements quantitatively confirm the overall validity of the 
effective potential of a 22-pole ion trap, which scales as r^°. For large rf amplitudes, 
however, new features in the potential have been observed in the form of ten almost 
equally spaced potential minima. These minima arise from the breaking of the 22- 
fold symmetry of the trap. They can be quantitatively explained by a slight tilt of 
half of the multipole rf electrodes within their mechanical tolerances. Also for other 
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multipole ion traps n — 1 minima in the effective trapping potential have been found 
as a consequence of a broken symmetry of the trap. 

This observation, which has become possible due to the high sensitivity of 
our photodetachment tomography scans, has implications for other spectroscopic 
experiments in 22-pole ion traps. In particular at cryogenic temperatures, trapped 
ions will reside predominantly in the ten pockets of the potential. Correspondingly, 
the ion density along the symmetry axis of the trap would become very small and 
only a small spectroscopic signal would be detected for a laser beam pointing along 
the trap axis. Further studies are needed to find out if ions that reside in the ten 
pockets may be subject to enhanced radiofrequency heating, similar to the influence 
of the static end electrode potential, which can increase the translational temperature 
by a few Kelvin [21]. Generally, it is therefore advisable to operate the ion trap at 
rf amplitudes far below 100 V to significantly suppress the pockets. To overcome the 
pockets a significantly enhanced precision in the manufacturing and assembly of 22- 
pole traps is required. An interesting alternative for precisely controllable multipole 
ion traps are planar, chip-based traps [24] , 
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Figure 1. (Color online) a) Measured density distribution of trapped OH~ 
ions at 300 K buffer gas temperature with the rf amplitude set to 160 V. The 
sketched geometry shows the layout of the ion trap, viewed along its symmetry 
axis. It includes the copper housing, the 22 rf electrodes (end-on), a surrounding 
shaping electrode, and the end electrodes, b) Zoom into the measured ion 
density distribution, each pixel represents an individual decay rate measurement, 
c) One-dimensional cut through the density distribution along the horizontal 
axis, d) Effective potential derived from the density distribution by assuming 
a Boltzmanm distribution of the trapped ions at 300 K. 
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Figure 2. (Color online) a) Ion density distribution at 170 K with an rf amplitude 
of 160 V. It shows a substructure of ten clearly distinct maxima, b) Cut through 
the effective potential, derived from the the density distribution. Overall, the 
potential is in accordance with the effective potential of an ideal 22-pole (solid 
line) . 




Figure 3. (Color online) a) Ion density distribution at 300 K with an rf- 
amplitude of 270 V. The histogram shows ten distinct maxima, which correspond 
to minima in the effective potential, b) The depth of the deepest minimum shows 
a strong increase as a function of the applied rf-amplitude. The solid line shows 
a polynomial fit with only a constant and a quadratic term. 
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Figure 4. (Color online) a) Depth of the deepest potential minimum as a function 
of the tilt angle of the rf electrodes with a potential applied to the end electrodes 
(upper points) and without any potential on the end electrodes (lower points). 
The solid lines show quadratic fits without a linear term. The inset shows the 
geometry of the 22-pole ion trap, viewed in the direction onto the coldhead (see 
Fig. [T^). The right wall of the trap together with the 11 implanted rf electrodes 
has been tilted by an angle of 1.0°. b) Calculated ion density distribution for 
300 K in the effective potential of 270 V rf amplitude and -2 V static end electrode 
potential. A tilt angle of 0.15° is chosen, which leads to a good agreement with 
the measured density distribution shown in Fig. |3^. 



